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Abstract 

We perform a Dalitz plot analysis of D + — > i"C~7r + 7r + decay with the CLEO-c data set of 
572 pb" 1 of e + e~ collisions accumulated at the ■0(3770). This corresponds to 1.6 x 10 6 D + D~ 
pairs from which we select 140793 candidate events with a small background of 1.1%. We compare 
our results with previous measurements using the isobar model. We modify the isobar model with 
an improved description of some of the contributing resonances, and get better agreement with 
our data. We also consider a quasi-model-independent approach and measure the magnitude and 
phase of the contributing Kir S wave in the range of invariant masses from the threshold to the 
maximum in this decay. This gives an improved descriptions of our data over the isobar model. 
Finally we allow for an isospin-two ir + ir + S wave contribution, and find that adding this to both 
the isobar model and the quasi-model-independent approach gives the best description of our data. 



2 



I. INTRODUCTION 



In comparison to other D + decay modes, the D + — > K n + n + decay is unique in many 
aspects. The large branching fraction, B(D + -> K-n + n+) = (9.51 ± 0.34)% [1], for this 
Cabibbo favored mode makes it the usual choice for normalization of other D + -meson decay 
rates. Understanding its peculiar intermediate substructure will be beneficial. The only 
obvious contribution to this decay, observed in the K7r-mass spectrum, is K*(892) n + , 
which comprises merely 12% of the total rate [1]. A large contribution of over 60% from a 
Kn S wave intermediate state has been observed in earlier experiments, including MARK 
III [2], NA14 [3], E691 [4], E687 [5], and E791 [6, 7], where the D + -> K-tx + tx + decay has 
been studied with the Dalitz plot technique [8]. Hence, the D + — > K~7t + tt + decay is a good 
laboratory to study Kn S wave dynamics. 

The previous analysis by E791 [6] achieved good agreement with their data by including 
a low- mass K~n + scalar resonance k that significantly redistributed all fit fractions (FF) 
observed by earlier experiments. This particular model, even though it is based on the 
largest data set, greatly disagrees with previous analyses and has been excluded from the 
average given by the Particle Data Group (PDG) [1]. 

There has been significant theoretical interest in this decay, sparked by the large, low- 
mass Kir S wave contribution. In Refs.[9-ll] the authors reanalyze the E791 data with 
their own models. E791 later reinterpreted their own data with a model-independent partial 
wave analysis [7], and we apply this in our analysis with minor modifications. 

The two identical pions in the final state should obey Bose symmetry. Assuming that the 
three-body decay is dominated by two-body intermediate states, there would be two identical 
K~7i + waves interfering with each other. This two-fold symmetry significantly reduces the 
degrees of freedom in the regular Dalitz plot analysis and allows the application of a model- 
independent partial wave analysis [7]. We would also expect a small contribution from 
the isospin-two S wave, which exhibits nontrivial dynamics as observed in scattering 

experiments [12]. 

The data used in this analysis were accumulated with the CLEO-c detector [13]. Our 
event sample is based on 572 pb" 1 of e + e~ collisions at y/s ~ 3774 MeV, produced by 
the Cornell Electron Storage Ring (CESR). This sample corresponds to the production of 
1.6 x 10 6 D + D~ pairs in the process e+e~ -> ^(3770) -> D+D-. We select 140793 D+ -> 
K~7t + tt + candidates for the Dalitz plot analysis (charge conjugation is implied throughout 
this paper). Our sample is very clean with a background fraction of about 1.1% and is 9 
times larger than the data set used by E791. The invariant mass resolution in this three- 
track D-meson decay is very good; we estimate it is better than 5 MeV/c 2 in most cases. It is 
improved by a kinematic fit requiring a three-track common vertex with the D-meson mass 
constraint. Our kinematic conditions are similar to those of MARK III, where D mesons 
are produced with small momentum. 

In Sec. II we briefly discuss CLEO-c experimental techniques, giving the event selection 
for the Dalitz plot analysis, the general fit method and methods to parametrize the signal 
efficiency and background distribution across the Dalitz plot. The formalism we use for the 
amplitude parametrization in this analysis is described in Sec. III. In Sec. IV we compare our 
results with the best previous measurements by E791 [6] and try to improve the isobar model 
in order to get a better description of our data. Finally, we apply a quasi-model-independent 
partial wave analysis, following Ref. [7] and measure the partial waves contributing to this 
decay in Sec. V. Systematic studies and cross-checks are considered in Sec. VI. We discuss 
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results and outstanding issues of this analysis in Sec. VII and summarize our results in 
Sec. VIII. In the appendix we discuss the kinematic variables and angular distributions 
used in this analysis. 



II. DETECTOR AND EXPERIMENTAL TECHNIQUE 

A. Detector 

CLEO-c is a general purpose detector which includes a tracking system for measuring 
momenta and specific ionization of charged particles, a ring imaging Cherenkov detector to 
aid particle identification, and a Csl calorimeter for detection of electromagnetic showers. 
These components are immersed in a magnetic field of 1 T, provided by a superconducting 
solenoid, and surrounded by a muon detector. The CLEO-c detector is described in detail 
elsewhere [13]. 

B. Event Reconstruction 

We reconstruct the D + — > K~7r + 7r + decay using three tracks measured in the tracking 
system. Charged tracks satisfy standard goodness of fit quality requirements [14]. Pion and 
kaon candidates are required to have specific ionization dE/dx in the main drift chamber 
within 4 standard deviations of the expected value at the measured momentum. 

In order to select D + — > K~it + it + decays, we use two kinematic variables, 



where E heara is the beam energy, and E D and P D are the energy and the magnitude, re- 
spectively, of the momentum of the reconstructed D + candidate. The tubc and \AE\ two- 
dimensional distribution and the projections for data are shown in Fig. 1. The resolutions 
in AE and m-Qc are represented as a(AE) = 6 MeV and a{m-Qc) = 1.5 MeV/c 2 , respec- 
tively; fits with a Gaussian function to the AE and mec peaks evaluate the resolutions to 
be 5.539 ± 0.014 MeV and 1.410 ± 0.013 MeV/c 2 , respectively. We require the events to fall 
in the "signal box" that is the overlap region of the AE and tubc signal regions defined as 
\AE\ < 2a(AE) and \m-Qc — ™>d\ < 2o"(mBc), respectively. In the case of multiple D-meson 
candidates per event we select the one with the smallest |Ai?| value. 

The K^n + n + final state has two identical n + mesons. To account for this symmetry we 
analyze events on the Dalitz plot by choosing x = m 2 (K~n + )i ow and y = m 2 (K~7r + ) high as 
the independent (x, y) variables. This choice folds all of the data onto the top half of the 
kinematically allowed region, as is shown in Fig. 2(a). The third variable z = m 2 (ii + 7i + ) 
is dependent on x and y through energy and momentum conservation. The invariant mass 
resolutions, propagated from the track error matrices, are shown in Fig. 3, and in 95% of 
cases are better than 5 MeV/c 2 . We use a kinematic fit to all 3-track candidates which 
enforces a common vertex and D + mass [1] constraint. We require that all events pass 
the kinematic fit successfully but do not restrict their x 2 . The kinematic-fit-corrected 4- 
momenta of all 3 particles are used to calculate invariant masses for further Dalitz plot 
analysis. Within its finite accuracy, the kinematic fit improves the K~it + -n + invariant mass 



AE = E D — -Ebeam, 



(1) 




(2) 
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resolution by 2 orders of magnitude. Proportional improvement is expected for all two-body 
invariant mass resolutions. 

After all requirements, we select 140793 events for the Dalitz plot analysis. The signal 
fraction in this sample / sig is estimated to be (98.917 ± 0.013)% from the fit to the mec 
distribution, shown in Fig. 1(b). In this fit the signal and background shapes are described 
by the double-Gaussian and ARGUS [15] functions, respectively, with all parameters free. 
This value of / s ; g is consistent with one obtained from the fit to AE distribution, shown in 
Fig. 1(c). In most fits to the Dalitz plot we use the fixed value of the signal fraction. Figure 2 
shows the Dalitz plot data and two projections onto the m 2 (K7r) axes [two entries per event 
for m 2 (K7r)\ ow and m 2 (-K7r)hi g h, respectively] and m 2 (im). Besides the clear if* (892) signal, 
no other narrow features are obvious. The strong left-right asymmetry of the if* (892) 
population density on the Dalitz plot is evidence of the interference between P and S 
waves. There are broad structures, including a peak at m 2 (Kn) around 1.3 GeV 2 /c 4 , a 
dip at m 2 (K7r) around 2.25 GeV 2 /c 4 , a dip at m 2 (7T7r) around 1 GeV 2 /c 4 , and a peak at 
m 2 (7T7r) around 1.6 GeV 2 /c 4 , that do not obviously correspond with known resonances or 
their reflections from other axes. These structures also do not correspond to a flat phase 
space distribution for nonresonant decays, since our efficiency is essentially flat across the 
Dalitz plot, and the background is very low. Thus, we are compelled to consider Kir strong 
interaction dynamics to explain the Dalitz plot. 
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FIG. 1: Event selection, (a) The mec and \AE\ two-dimensional distribution for data and the 
projections on (b) ttibc and (c) AE. The mec and AE signal regions, defined in the text, are 
shown as the bands in the figures. In (a), the "signal box" is indicated as the crossing area of 
the two bands while the "sideband box," defined in the text, is indicated as the shaded rectangle. 
Each projection is made with the events in the signal region of the other kinematic variable; the 
fit curve, described in the text, is also shown. 
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FIG. 2: (a) Dalitz plot for data and their projections on (b) m 2 (K7r) (two entries per event), and 
(c) m 2 (7T7r) variables. 
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FIG. 3: Invariant mass resolutions before the kinematic fit for (a) m(KTr)hi g h, (b) m(Ki:)i ow , and 
(c) m(7T7r). 

C. Fit method 

In order to describe the event density distribution on the Dalitz plot we use a probability 
density function (p.d.f.), V(x,y), which depends on the event sample being fit: 

!J\f e e(x,y) for efficiency, 

NbB{x,u) for background, 

f sig Af s \M{x,y)\ 2 e(x,y) + (1 - f sig )Af B B(x,y) for signal with background, 

where the e(x,y) and B(x,y) are the functions representing the shape of the efficiency and 
background, respectively, across the Dalitz plot. The signal p.d.f. is proportional to the 
efficiency-corrected matrix element squared \Ai(x,y)\ 2 , defined in Sec. Ill, whose fraction 
/si g is introduced earlier. The background term has a relative (1 — / S j g ) fraction. All p.d.f. 
components are normalized separately using the normalization integrals over the Dalitz plot 
area l/J\f E = J e(x,y)dxdy, 1/Mb = J B(x,y)dxdy, and 1/Ms = J \M{x,y)\ 2 e(x,y)dxdy, 
which provides the overall p.d.f. normalization, / V(x,y)dxdy = 1. The p.d.f. free parame- 
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ters are optimized with an maximum likelihood fit that minimizes the sum over N events: 

N 

C = -2j2^ogV(x n ,y n ). (4) 

n=l 

To estimate the fit quality we use Pearson's statistics for adaptive bins, similar to our 
previous analysis [16]. 

D. Efficiency Parametrization 

To determine the efficiency we use a signal Monte Carlo (MC) [17] simulation where one of 
the charged D mesons decays in the signal mode uniformly in phase space, while the other D 
meson decays in all known modes with relevant branching fractions. These underlying events 
are input to the CLEO-c detector simulation and processed with the regular reconstruction 
package. The MC-generated events are required to pass the same selection requirements as 
data selected in the signal box, as shown in Fig. 1(a). In each event we consider only the 
signal mode side to prevent nonuniformity of the efficiency due to the resonance substructure 
of the other side D decay currently implemented in our generic simulation. The efficiency of 
the Ktiti final state selection for the Dalitz plot analysis is estimated to be (51.11 ± 0.07)% 
where the error is only statistical. This number also accounts for a correction factor, 0.984, 
due to the nonuniform population of the data on the Dalitz plot, 

j = gData = En=l e ( X n,yn) / N 

e M c I e(x,y)dxdy / J dxdy 

where £oata and £mc are the average efficiencies for data and uniformly generated MC sam- 
ples, the function e(x,y) is an efficiency over Dalitz plot defined later by Eqs. (6) and (7), 
the sum runs over all N events in the data sample, and the integrals are taken over the area 
of the Dalitz plot. 

To parametrize the efficiency e(x, y), we use a symmetric third-order polynomial function 
centered on the arbitrary point (x c , y c ) = (1.5, 1.5) (GeV/c 2 ) 2 on the Dalitz plot. With 
x = x — x c and y = y — y c , the efficiency is the product of the polynomial function: 

e(x, y) = T(v)[l + E 1 (x + y) + E 2 (x 2 + f) + E 3 (x 3 + f) + E xy xy + E xyn (x 2 y + xy% (6) 

and sinelike threshold factors for each Dalitz plot variable v (= x,y or z): 

T ( v \ = J sin(£th,t;k - Umaxl), for < E thjV \v - v mSLX \ < tt/2, , 7 
1 ' \ 1, for E thtV \v - v max \ > n/2, 1 ' 

where all polynomial coefficients Ei, E 2 , E 3 , E xy , E xyn , and E t ^ v are the fit parameters. 
Each variable v has two thresholds v m i n and v max . We expect low efficiency in the regions 
v v ma x only, where one of three particles is produced with zero momentum in the D-meson 
rest frame and thus has a small momentum in the laboratory frame. Figure 4 and Table I 
show results of the fit to the entire signal MC sample of D + — > K~ir + ir + events selected on 
the Dalitz plot. The polynomial function with threshold factors well describes the efficiency 
shape. If we consider subsamples of our signal MC, such as D + versus D~ , we find that the 
variation of the efficiency parameters is small compared to their statistical uncertainties. In 
fits to data we use this efficiency shape with fixed parameters, and variations constrained 
by the errors from our fit to the signal MC are allowed as systematic checks. 
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TABLE I: Fit parameters for the efficiency from the signal MC sample. 



Parameter 


Value 




-0.0153±0.0090 




-0.030±0.011 


E3 


0.f62±0.020 


E 


-0.053±0.019 


E 


0.673±0.055 


Eth,x = Eth,y 


4.25±0.23 


Eth,z 


2.907±0.075 


Pearson x 2 l v 


649/573 


Probability (%) 


1.5 


Events on DP 


477978 




FIG. 4: For the efficiency shape: (a) Dalitz plot of the signal MC generated uniformly in phase 
space and its projections on (b) m 2 {Kir) (two entries per event) and (c) m 2 (7T7r) variables. The solid 
histogram is a projection of the function described in the text which parametrizes the efficiency. 
Binned results are shown, but the efficiency shape is determined with an unbinned maximum 
likelihood fit. 



E. Background Parametrization 

A shape for the background on the Dalitz plot is estimated using data events from a mec 
sideband region, shown by the hatched box in Fig. 1(a). This box is shifted in AE from the 
signal region to have the same K~tt + tt + invariant mass range as candidates in the signal box. 
We consider only events from the low-mass mec sideband as the high-mass sideband has a 
significant contribution from signal events due to a "tail" caused by initial state radiation. 
This tail is clearly seen in the mec distribution shown in Fig. 1(b). 

The background is a small contribution that has little effect on our fits. Nevertheless, 
we study the background composition using generic MC simulation for all known modes 
and find the following. The pileup of events at m 2 (K 71)^1^ ~ 2.6 GeV 2 /c 4 is caused by the 
misreconstructed D decays from D°D° and D + D~ pairs and the combinatorial background 
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from the e + e~ — > gg (q = u, d and s) continuum; their contributions have the relative 
fractions of 62%, 13% and 25%, respectively, in this area. The dominant misreconstructed 
D decays are D° — > K~-n + it°, ai(1260) + K~ , and D + — > K~it + it + , where one of the final 
state pions is misreconstructed and replaced with a n + meson from the other D decay. The 
shape of the background is well reproduced in our simulation for events selected from either 
the signal or sideband box. 

To parametrize the background shape on the Dalitz plot we employ a function similar to 
that used for the efficiency, shown in Eqs. (6) and (7). Figure 5 and Table II show results 
of the fit with the background polynomial function to our sideband sample. In cross-checks 
with subsamples we find the variation of the shape parameters is small compared to their 
statistical uncertainties. We use the nominal background shape with fixed parameters in fits 
to data, and allow the parameters to vary constrained by their errors as a systematic check. 

TABLE II: Fit parameters for the background shape from the fit to the sideband region. 



Parameter 


Nominal value 


Bi 


0.63±0.22 


B 2 


0.95±0.39 


B 3 


0.41±0.54 


B X y 


-0.20±0.62 


Bxyn 


-1.2±1.3 




1.31±0.13 


Bth,z 


11.2±6.5 


Pearson x 2 /v 


129 / 97 


Probability (%) 


1.6 


Events on DP 


1554 



III. DECAY AMPLITUDE PARAMETRIZATION 
A. Matrix element 

In this analysis we follow the formalism of E791 [7] with only minor variations. In 
the formulas below for the Dalitz plot variables, we also use Mandelstam notations s = 
m 2 (K^7Ti ), t = m 2 (fT _ 7T2~), and u = m 2 (7r ] f n})- We choose two of them s and t as 
independent, and the third u is dependent, constrained by the equation s+t+u = m 2 D +m 2 K + 
2m 2 . Then, the matrix element has an explicit Bose-symmetric form for pion permutations 

M(s, t) = A(s, t) + A(t, s) + A[= 2 ( M ( S , t)). (8) 

Below we discuss the amplitudes contributing to the matrix element. 
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FIG. 5: (a) Dalitz plot of data in the sideband box and projections on (b) m 2 (K7r) (two entries 
per event), and (c) m 2 (7T7r) variables. The solid histogram shows the projection of the fit function 
used to parametrize the background shape described in the text. 
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FIG. 6: Three-body decay d — > Rc — > abc in the resonance R rest frame. 



B. Partial Ktv amplitudes 



Each Ktt amplitude is defined using a sum over the decay orbital momentum L of two- 
body partial waves 

A(s, i)= £ n L (s, t)rt(q(s))A L (s), (9) 

L=0 

with parameters as described below. In this analysis we consider the sum up to the maximal 
orbital momentum L max = 2. 

We assume the D + — > K~7t + tt + decay goes via a quasi-two-body intermediate state, 
d — > i?c, containing the resonance i? and particle c, followed by the decay of the resonance 
to the final stable particles a and b, R — > a&. This is shown schematically in Fig. 6. The 
term in Eq. (9) represents the angular distribution, which we use in the invariant 
forms [18] 

n L=0 (m 2 ab ,m 2 ac ) = 1, (10) 
fiL=i(m afe , m ac ) = m 5c - m ac + , (11) 



10 



V L =2(ml b , ml) = [n L=l f - \ (ml - 2m\ - 2m\ + K m c? 



" V m ab 

2 fm? — m 2 ^ 2 



ximl- 2mi - 2ml + ^7^" ' (12) 

V m ab ' 

where m^, m a , mj, and m c are the masses of decaying and product particles and m a i>, m ac , 
and mbc are the relevant invariant masses. In the appendix we show that these angular dis- 
tributions are equivalent to those applied in the E791 analysis [7] up to constant coefficients. 

The form factors T^{q) in Eq. (9) and ^(q) in Eqs. (20) and (22) are defined using the 
Blatt-Weisskopf form [19] 

L = 0: JF°(<?) = 1, (13) 
L = l: ^(q) 



(14) 



L = 2: ^(q) = 



W + d (15) 



\ 9 + 3g 2 + g 4 



where the index V stands for the D or R decay vertex, q = r v P, P is the magnitude of 
the momentum of the decay products in the decaying particle's rest frame, and ry is the 
effective radius for the D or R vertex, respectively. For both D and R decays, q v = r v P v , 
where Py is the magnitude of momentum of the decay products calculated at m ab = m R , 
the pole mass of R. The form factors are normalized by the condition J-y{qv) = 1- 

The values of radial parameters are discussed in Sec. IV. Expressions for the decay 
products' momentum for both vertices can be found in the appendix. The mass dependences 
of J~p(q(s)) form factors for resonances are shown in Table III. 

Another Gaussian form factor shape of the scalar resonance, 

J$(q) = e-ta 2 -*)/" (16) 

is applied in the E791 [7] analysis. This shape is suggested by Tornquist [20], and has a steep 
dependence on momentum. A reanalysis of the E791 data [10] found that this form factor is 
not required by the data. We use only the Gaussian form factors J-^ and JF^ from Eq. (16) 
for the scalar components k and 1^(1430), when comparing results with E791 model C [6]. 
We use the unit form factor from Eq. (13) for scalar resonances in all other models and with 
the binned partial waves discussed below. 

The partial waves Al(s) in Eq. (9) are the angular momentum L-dependent functions 
of a single variable z, which is either s or t. In the D + — > K~it + it + decay the S, P and 
D waves (L = 0, 1, or 2, respectively) are represented by the sum of functions Wr for 
individual intermediate states 

+ W S ,binned; 

(17) 

Al(z) = Wk*(S92) + Wk* (1410) + Wk* (1680) + Wp.binned, (18) 

A 2 (z) = W K 

|(1430) 

+ w D 

, binned' 

(19) 

The contribution of nonresonant decays is represented by c^r = aNRe l ^ NR , a complex con- 
stant with two fit parameters for magnitude clnr and phase 4>nr- The Wx,,binned are the 
binned amplitudes as discussed below. 

W R = CrWr^rP) (20) 
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is the shape of an individual resonance Wr [see Eqs. (21)-(24)] multiplied by the form 
factor in the resonance R decay vertex J~r{trP) and the coupling constant cr = aRe l ^ R . 
The resonance R production magnitudes cir and phases 4>r are parameters of the fit to the 
Dalitz plot. 



C. Resonance shapes 

For intermediate Kn resonances we use the standard Breit-Wigner function 

yV R (m) = — — -, (21) 

m R — m z — trriRL (m) 

where m 2 = s and the mass-dependent width has the usual form 

r(m) = r^(-^) 2L+1 [^(r R P)] 2 . (22) 

For _K"g(1430) we have tested both the Breit-Wigner function [Eq. (21)] and the Flatte 
parametrization 

Wfl(m) = — 2 4^-2 7 — \ > ( 23 ) 

m R -m 2 -l Eab 9RabPab{m) 

where qr^ is a coupling constant of resonance R to the final state ab and p a b(m) = 2P/m is 
a phase space factor. We test the Flatte parametrization because the i^o(1430) mass is close 
to the Krj and Krj' thresholds, which could significantly distort the resonance shape [10]. 
and the complex pole proposed in Ref. [11] 

W R (m) = l —, (24) 

s R — m z 

where sr is a pole position in the complex s = m 2 plane. This function represents the 
first term of the Laurent series in the expansion of the chiral perturbation theory complex 
amplitude for the scalar wave. This approach is common, and a survey of pole positions 
extracted from different experiments can also be found in Ref. [9]. This complex pole is 
equivalent to a Breit-Wigner function with constant width. 



D. Isospin-two 7T+7T+ S wave amplitude 

The isospin-two S wave amplitude in Eq. (8) is a sum of two components 

A^Ku) = c^Wt 2 + W^Led, (25) 

where c W7r = a^e' 1 ^ is a complex coupling constant and Wgjj? nned is discussed in Sec. HIE. 
The first term of this sum is parametrized by a unitary form [21] 

Wf V) = " g( '" )e !" m) " -, (26) 

where m is a invariant mass, rjl(m) is an inelasticity, and S^m) is a phase of the n + n + 

wave with total spin and isospin two. The phase <5o( m ) is assumed to be proportional to 
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the decay momentum at threshold and sculpted by a polynomial function at higher mass 
range 



W = - 



—a^Jm 2 /4 — m 2 



bm? 



cvrr 



(27) 



with parameters a = (55.21 ± 3.18) deg/GeV, b = (0.853 ± 0.254) GeV" 2 , c = (-0.959 ± 
0.247) GeV~ 4 , and d = (0.314 ±0.070) GeV~ 6 , obtained in Ref. [21] from a fit to the data of 
the scattering experiments [12]. We use this function with fixed parameters. The inelasticity 
r] 2 (m) in the mass range of m < m min ps m(p + p + ) pal. 5 GeV/c 2 is expected to be near unity. 
Then f] 2 (m) decreases due to the 7t + tt + — > p + p + rescattering at a higher mass range. In our 
fits we use a smooth approximation for this threshold behavior 



with fit parameters m r 



r i, 
i 
i 



2 



COS 7T 



m, 



iax , and Arj. 



m max m min 



m < m min 
, m min <m<m B 
m > m max , 



(28) 



E. Binned amplitude 

The complex term W^bmned m Eqs. (17)— (19), and (25), where L=0, 1, or 2, is intended 
to provide a completely model-free parametrization of the partial wave. It can be used alone 
or in combination with other terms. In the latter case it represents a correction to the 
complex amplitude of the isobar model. We use this term in the form of an s-dependent 
complex number 

W^bin„cd(s) =a L (sy Ms \ (29) 

with functions ol(s) and 0l(s) defined by an interpolation between the bins for the mag- 
nitude, ciLk, and phase, (pLk, where k(s) = 1,2, ...,Nl is an s-dependent index of these 
bins. For all K7r waves we define Nl = 26 uniform bins in s = m 2 Kn in the range [0.4,3.0] 
(GeV/c 2 ) 2 , as shown in Table III. Similar, for I = 2 7T+7T+ S wave we define iY"f= 2 = 18 
uniform bins in u = [s in Eq. (29) is replaced with u] in the range [0.1,1.9] (GeV/c 2 ) 2 . 
This binning scheme covers the kinematically allowed range of the Kir [0.633,1.730] GeV/c 2 
and 7r + 7T + [0.279,1.376] GeV/c 2 invariant mass spectrum in the D + — > K~it + it + decay. We 
interpolate linearly between bin centers in our fitting function. 



F. Fit fraction 

We estimate a contribution of each component in the matrix element using a standard 
definition of the fit fraction, 

FF J\A R (x,y)\ 2 dxdy 

R f\M(x,y)\ 2 dxdy> 1 ] 

where Ar(x, y) is an amplitude contribution from the R component to the total matrix 
element M(x,y) from Eq. (8) and the integrals are taken over the area of the Dalitz plot. 
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TABLE III: Bins for the Ktt S, P, and D waves and Blatt-Weisskopf form factors for if* (892) and 
K*(1680) from Eq. (14), #2(1430) from Ec l- ( 15 )> and ^o( 1430 ) from Eq. (16), calculated in the 
D-meson decay vertex. 



Bin m 


i (GeV/c 2 ) 2 

Kit \^"^ " 1 ° 7 


m Kn (GeV/c 2 ) 


Blatt-W. form factors J-^(q) for 


Gaussian FF 


# 

7T | 


Bin range 


Bin range 


U enter 


Z^* / or^O^ 

A (092J 


K ^IboUJ 


is* l 1 a on\ 
A 2 (1430) 


Krt(1430) 


1 


0.4 — o.5 


n coo n TriT 

0.632 — 0.707 


n C71 
U.DYl 


0.888 


n OFH 

U.zoU 


n one 

U.oUo 


0.347 


2 


0.5 — 0.6 


0.707 — 0.775 


n T A o 

U.742 


n m o 

0.918 


n o F;n 

U.zo9 


n oo /i 

U.oz4 


0.380 


3 


0.6 — 0.7 


0.775 — 0.837 


0.806 


n n /i o 

0.948 


n oeT' 

U.zbY 


U.o4o 


0.415 


4 


0.7 — 0.8 


n o o t r\ on a 

0.837 — 0.894 


0.866 


n noo 

0.982 


n ott 

0.277 


n oco 

0.368 


0.451 


5 


0.8 — 0.9 


n on a r\ n -in 

U.o94 — U.94y 


n noo 

U.9zz 


1.U17 


n oot 

U.zo7 


n on /I 

U.o94 


0.489 


6 


0.9 — 1.0 


n n/in i nnn 

u.y4y — i.uuu 


n nTF; 

u.y (o 


l.Uoo 


n oo7 
U.zy 7 


n a oi 
U.4zl 


0.528 


7 


1.0 — l.l 


1 nnn i n/in 
l.UUU — l.U4y 


1 no 
l.Uzo 


l.Uyo 


n onn 

u.ouy 


n /i f;o 
U.4oz 


0.570 


8 


1.1 — 1.2 


1 n/in i nnr; 

l.U4y — l.Uyo 


1 nTo 
l.U/z 


i i /i n 
1.14U 


n 001 
U.ozl 


n /I of; 

U.4oo 


0.612 


9 


1.2 — 1.3 


1 nn k i i /i n 
l.Uyo — 1.14U 


1 1 1 Q 

l.llo 


1 1 oo 
l.loo 


n oof; 
U.OOO 


n f;oo 
U.OZO 


0.656 


10 


1.3 — 1.4 


1.14U — l.loo 


1 1 £?0 

l.lbz 


i o a n 

l.z4U 


n o /i n 

0.349 


n ft c a 

0.564 


0.700 


11 


14 — 1.5 


1 1 OO 1 oo fc 

l.loo — l.ZZO 


1 on a 
l.zU4 


1 one 

1.29b 


n oc n 

0.365 


0.609 


0.746 


12 


1.5 — 1.6 


1 oof; i o/^f; 
l.zzo — l.ZOO 


i o/t f; 
1.z4o 


1 o f;o 
l.OOO 


0.3o3 


0.659 


0.792 


13 


1.6 — 1.7 


l.ZOO — l.oU4 


1 oof; 

l.zoo 


1 A o K 
1.4Z0 


n /i no 
0.4U2 


0. ( 15 


0.839 


14 


1.7 — is 


1 on A 1 O A o 
l.oU4 — - l.o4z 


1 ooo 
l.OZO 


1 /i nn 
1.4yy 


n /loo 
0.423 


n "tip. 
0. ( /6 


0.885 


15 


1.8 — 1.9 


1 O /l O 1 O TO 

1.342 — 1.378 


1.360 


1 CT01 

1.581 


ft A An 

0.446 


n o /I /i 

0.844 


0.932 


16 


1.9 — 2.0 


1 ^78 1/11/1 
l.O/o 1.4fcl4fc 


1 QQR 

i.oyo 


1 P.79 
1.0 1 Z 


n /1 71 


n Qis 
u.y io 


0.978 


17 


2.0 — 2.1 


1.414 — 1.449 


1.432 


1.773 


0.500 


0.999 


1.023 


18 


2.1 — 2.2 


1.449 — 1.483 




x.oou 




l.UOu 


1.067 


19 


2.2 — 2.3 


1.483 — 1.517 


1.500 


2.013 


0.568 


1.178 


1.109 


20 


2.3 — 2.4 


1.517 — 1.549 


1.533 


2.157 


0.608 


1.275 


1.150 


21 


2.4 — 2.5 


1.549 — 1.581 


1.565 


2.320 


0.654 


1.374 


1.189 


22 


2.5 — 2.6 


1.581 — 1.612 


1.597 


2.506 


0.707 


1.472 


1.226 


23 


2.6 — 2.7 


1.612 — 1.643 


1.628 


2.719 


0.766 


1.568 


1.259 


24 


2.7 — 2.8 


1.643 — 1.673 


1.658 


2.962 


0.835 


1.657 


1.290 


25 


2.8 — 2.9 


1.673 — 1.703 


1.688 


3.240 


0.913 


1.737 


1.318 


26 


2.9 — 3.0 


1.703 — 1.732 


1.718 


3.554 


1.002 


1.805 


1.342 


rriK-K, 


min = %- + = 0.633 G 


reV/c 2 










rriKn, 


max — m D+ ~ 


rjVH = 1.730 C 


r eV/c 2 











G. Expected contributions 

A priori, all known Kn resonances in the mass range from the production threshold 
up to 1.73 GeV/c 2 , such as if* (892), fT*(1410), AT* (1430), fT*(1430), AT*(1680), and even 
higher mass resonances, might contribute to the amplitude of the D + — > K~n + n + decay. 
Table IV shows parameters of Kn resonances which have been considered in this analysis. 
One would also foresee an / = 2 n + n + S wave final state interaction. Simulations of some 
of the expected contributions to the Dalitz plot are shown in Fig. 7. In contrast to data 
[Fig. 2(a)] the K*(892) population density in Fig. 7(b) is symmetric without interference 
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with an S wave. 



3540108-007 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 



m 2 (K7i) low (GeV/c 2 ) 2 

FIG. 7: Simulation of the expected contribution to the D + — ► iT~7r + 7r + Dalitz plot from various 
intermediate states, (a) Low mass Kir S wave (k), (b) K*(892)tt, (c) K(J(1430)tt, (d) iq(1430)vr, 
(e) ^*(1680)7r, and (f) K{irit) I=2 with 1 = 2 7r + vr+ 5 wave. 



IV. FITS TO DATA USING THE ISOBAR MODEL 
A. Comparison with E791 Model C 

First we compare our results, obtained in the framework of the isobar model, with E791 
models A, B, and C from Ref. [6]. In particular, the most complete model C contains 
iT*(892)7r+, K*(U30)n + , K*(U30)n + , K*(1680)n + , ktt + , and nonresonant (NR) contribu- 
tions. Following E791 we allow a scalar Kir amplitude, the " Breit- Wiener resonance 
with mass-dependent width. We set ck*(892) = 1 in Eqs. (18) and (20), and all other magni- 
tudes and phases are defined with respect to if* (892). Gaussian form factors ^,(ri)P c ) and 
jF^(r^P) from Eq. (16) are used for ^(1430) and k. For all Kir resonances with nonzero 
spin, the radii in the Blatt-Weisskopf [19] form factors ro = 5 GeV^ 1 and tr = 1.5 GeV -1 
are fixed to the values used by E791. On the Dalitz plot, the p.d.f. for this model looks 
indistinguishable from statistics shown in Fig. 2(a). The Dalitz plot projections with p.d.f. 
components are shown in Fig. 8. Fit fractions and phases, obtained in our fit, are compared 
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with E791 [6] in Table V and are statistically consistent. Magnitudes are not comparable 
because of a different choice for normalization. The values obtained for resonance parame- 
ters are compared in Table IV. In particular, we get dominant contributions from S wave 
components: the NR, Kq(1430)it, and kit fit fractions are ss9%, 10%, and 33%, respectively. 
The total sum of all fit fractions is 65.5%, indicating substantial constructive interference. 
Apparently, model C gives a poor fit quality x 2 / z/= 53 1/391. A large discrepancy between 
the fit and the data is seen in Fig. 8 for the m 2 (7r + 7r + ) projection in the range of [1.4, 1.9] 
(GeV/c 2 ) 2 . That motivates us to explore alternative models of the decay amplitude. 

TABLE IV: The Breit-Wigner resonance parameters used or measured in the isobar model; CLEO- 
c vs E791. The fixed parameters used in the E791 experiment are taken from PDG 2000 [22] 
with their uncertainties shown in square brackets. The measured values are shown with two 
uncertainties: statistical and systematic. The values shown in parentheses (with statistical error 
only) were obtained in cross-checks when these parameters are allowed to float. 



Parameter 


E791 [PDG 2000] 


CLEO-c 


PDG 2006 [1] 


(MeV/c 2 ) 


Model C (if float) 


Model 12 (if float) 




m K*(892) 


896.1 [±0.27] 


896 (894.8±0.5) 


895.7±0.2±0.3 


896.00±0.25 


IV* (892) 


50.7 [±0.6] 


50.3 (45.5±0.4) 


45.3±0.5±0.6 


50.3±0.6 


"I** (1430) 


1459±7±12 


1463.0±0.7±2.4 


1466.6±0.7±3.4 


1414±6 


IV* (1430) 


175±12±12 


163.8±2.7±3.1 


174.2±1.9±3.2 


290±21 


m K*(1430) 


1432.4 [±1.3] 


1432.4 (1436±11) 


1432.4 (1427±7) 


1432.4±1.3 


rX*(1430) 


109 [±5] 


109 (132±21) 


109 (120±13) 


109±5 


m K*(1680) 


1717 [±27] 


1717 (1782±41) 


1717 (1679±59) 


1717±27 


IV*(1680) 


322 [±110] 


322 (565±131) 


322 (446±119) 


322±110 


m K*(U10) 


1414[±15] 


1414 


1414 


1414±15 


I 1 A"* (1410) 


232[±21] 


232 


232 


232±21 


m K 


797±19±43 


809±1±13 


Complex pole, 


i-C *(800) is not 


r K 


410±43±87 


470±9±15 


see Table VI 


established 



B. Variations of Model C 

The Gaussian form factors, given by Eq. (16) for scalar resonances, J r< p i (r R P), and 
jF^(r D P c ), behave similarly to the L = 4 Blatt-Weisskopf form factor. This behavior is 
not preferred by either our or the E791 data; see Ref. [10] for details. With our data we 
find that results are not very dependent on the assumed resonance decay vertex form factor 
jF^(r^P). However, the D-meson decay vertex form factor jF^r^P) changes the S wave 
dependence on s significantly. In particular, this Gaussian form factor (see Table III) sup- 
presses the contribution of Kit at low mass. To agree with the data when fitting with this 
factorized form factor, the magnitude of the complex function (S wave) increases at low 
Kir mass which gives an illusion of resonance behavior. For all models other than model C 
below we use unit form factors from Eq. (13) for S wave contributions. 

For the ^(1430) resonance we measure m K *( 1430 ) = 1463.0 ± 0.7 ± 2.4 MeV/c 2 and 
fV *(i430) = 163.8 ±2.7±3.1 MeV/c 2 , which are consistent with E791 results, but inconsistent 
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TABLE V: Comparison of CLEO-c results with E791 using the isobar fit, model C. Shown are the 
fitted magnitudes, a in arbitrary units, the phases, (f> in degrees, defined relative to the K*(892)ir + 
amplitude, and the FF. Also indicated are the fitted masses m and widths V of the spin-zero 
resonances. Magnitudes a are not comparable between the two experiments because of a different 
choice for normalization. 



Mode 


Parameter 


E791 


CLEO-c 


NR 


a 


1.03±0.30±0.16 


7.4±0.1±0.6 




<M°) 


-11±14±8 


-18.4±0.5±8.0 




FF (%) 


13.0±5.8±4.4 


8.9±0.3±1.4 


iT (892)vr+ 


a 


1 (fixed) 


1 (fixed) 




<M°) 


(fixed) 


(fixed) 




FF (%) 


12.3±1.0±0.9 


11.2±0.2±2.0 


Ko(1430)tt 4 


a 


1.01±0.10±0.08 


3.00±0.06±0.14 




<t>n 


48±7±10 


49.7±0.5±2.9 




FF (%) 


12.5±1.4±0.5 


10.4±0.6±0.5 




m (MeV/c 2 ) 
T (MeV/c 2 ) 


1459±7±12 


1463.0±0.7±2.4 




175±12±12 


163.8±2.7±3.1 


l?2(1430)7r 4 


a 


0.20±0.05±0.04 0.962±0.026±0.050 




<M°) 


-54±8±7 


-29.9±2.5±2.8 




FF (%) 


0.5±0.1±0.2 


0.38±0.02±0.03 


lT (1680)vr 4 


a 


0.45±0.16±0.02 


6.5±0.1±1.5 




<M°) 


28±13±15 


29.0±0.7±4.6 




FF (%) 


2.5±0.7±0.3 


1.28±0.04±0.28 




a 


1.97±0.35±0.11 


5.01±0.04±0.27 




<M°) 


-173±8±18 


-163.7±0.4±5.8 




FF (%) 


47.8±12.1±5.3 


33.2±0.4±2.4 




m (MeV/c 2 ) 
T (MeV/c 2 ) 


797±19±43 


809±1±13 




410±43±87 


470±9±15 


Form factor 


^(GeV- 1 ) 


1.6±1.3 


1.5(fixed) 




^(GeV- 1 ) 
isTTT ^(GeV" 1 ) 


5.0±0.5 


5 (fixed) 


Other R -► 


1.5 (fixed) 


1.5 (fixed) 


£ FF (%) 




88.6 


65.5 


Goodness 


X 2 /" 


46/63 


531/391 



with current PDG [1] values, as demonstrated in Table IV. Similar behavior is reported in 
Ref. [23] from the FOCUS Collaboration. In Ref. [10] Bugg surmises that the if* (1430) 
resonance parameters might change due to the opening of the Krj' channel. In order to 
accommodate this effect, we test the Flatte parametrization of Eq. (23), which depends on a 
floating mass mx*(i43o) and three coupling constants, Qk-ki 9kt), and Qk^'- We find that our 
data are consistent with g Kr) = 0, and this coupling is dropped from further consideration. 
The resulting values of the other parameters are shown in Table VI. We do not find any 
significant difference between Breit-Wigner and Flatte parametrizations in the shape of the 
Kq(1A30) complex amplitude or in the fit quality. 
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TABLE VI: Alternative parameters obtained in the fits with the isobar model. 



Mode 


Amplitude Parameter (MeV/c 2 ) 


Model C 


Model 12 


Kg(l430)7T 4 


' Breit-Wigner m K * {U30) 


1463.0±0.7±2.4 


1466.6±0.7±3.4 




F if* (1430) 


163.8±2.7±3.1 


174.2±1.9±3.2 




Flatte "ix*(i430) 


1462.5±3.9 


1471.2±0.8 




9Kn 


532.9±8.5 


546.8±4.2 




9K v 










9Kri' 


197±106 


230±32 




Breit-Wigner m K 


809±1±13 


888±2 




r K 


470±9±15 


550±12 




Complex pole $lm K 


769.9±6.3 


706.0±1.8±22.8 






-221.2±8.4 


-319.4±2.2±20.2 



In model C we also measure the k resonance Breit-Wigner parameters m K = 809 ± 1 ± 
13 MeV/c 2 and T K = 470 ± 9 ± 15 MeV/c 2 . In Ref. [11] Oiler emphasizes that the Breit- 
Wigner resonance with a mass-dependent width is not the best choice for the Kn low-mass 
phase parametrization. Following his prescription, we replace the Breit-Wigner function by 
the complex pole from Eq. (24) with the initial pole position at s K = (710 — i310) 2 MeV 2 /c 4 . 
We allow the pole position to float and obtain its optimal location, as shown in Table VI. 

We also test the JT*(892) parameters, as shown in Table IV. The K*(892) mass is 
consistent with the PDG [1] value, while the width is about 5 MeV/c 2 smaller. In further 
analysis we allow the if* (892) mass and width to float. 

Inclusion of a if* (1410) component in the fit does not result in any significant contribu- 
tion. 

C. Model 12 

Tuning of the isobar model for Kn waves only does not improve significantly the probabil- 
ity of consistency between the data and the model, which is still small. The large discrepancy 
in the m 2 (n + n + ) spectrum persists even if we use the model-independent parametrization 
for the Kn S wave considered below. We solve this problem by including in the matrix 
element [Eq. (8)] a contribution from the 1 = 2 n + n + S wave, using Eq. (26), which leads to 
a model called model 12. The threshold of the 7r + 7r + — > p + p + rescattering process is located 
at the upper edge of the m 2 (7r + 7r + ) kinematic limit in our decay. The edge of the kinematic 
border does not allow us to resolve correlations between the Arj and m max parameters in 
Eq. (28). We fix the value of Arj = 1, while floating the m min and m max parameters. In 
model 12, compared to model C, the unit form factors F% = jFjj, = 1 are used for the scalar 
components, the complex pole amplitude from Eq. (24) is used for the low-mass Kn S wave, 
and the K*(892) parameters are allowed to vary in the fit; other conditions stay the same 
as in model C. 

Model 12 gives the best analytical description of our data. The obtained fit parameters 
are presented in Table VII. The total p.d.f. and separate components are shown in the Dalitz 
plot projections in Fig. 9. The large discrepancy in the m 2 (n + 7i + ) spectrum is eliminated, 
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improving the probability of consistency between the model and data to 13%. 
D. Variations of Model 12 

We did not find any significant contribution from i^*(1410) by including it in the fit. 
This resonance is excluded from further consideration. 

We also test the alternative descriptions for the 1^(1430) parameters as shown Table VI. 
We do not find any significant difference between Breit-Wigner and Flatte parametrizations 
in the shape of the Kq(1430) complex amplitude or in the fit quality. We therefore use the 
Breit-Wigner function in model 12 with the mass and width of -K"q(1430) allowed to float. 

The test of alternative descriptions for the low-mass Kn S wave is presented in Tables VI 
and VII. The pole amplitude slightly improves the fit quality (Ax 2 = —10) compared to the 
Breit-Wigner resonance with mass-dependent width. This substitution leads to a significant 
redistribution of the NR and n components of the scalar wave, though their sum results in 
very small variation of the complex function. In model 12 we use the pole amplitude for the 
Kir low-mass S wave. 

V. QUASI-MODEL-INDEPENDENT PARTIAL WAVE ANALYSIS 
A. QMIPWA for Kir S wave 

The biggest issue of any Dalitz plot analysis is its model dependence. An attempt to 
mitigate the model dependence for the decay under study is described in [7]. Here we 
reproduce this analysis using a slightly modified technique, which we call the quasi-model- 
independent partial wave analysis (QMIPWA). We apply this technique as an extension of 
our model 12. 

In QMIPWA we modify the parameters of model 12 as follows. The NR and ktt compo- 
nents in model 12 are replaced by the binned S wave amplitude [Eq. (29)]. The 26 binned 
S wave magnitudes and phases of m 2 Kn are floating fit parameters. We keep the Kq(1430)tt 
contribution in its Breit-Wigner form, because it has a sharp structure, that cannot be well 
reproduced by a binned amplitude. The Kq(1430) parameters are fixed to their values from 
model 12 in order to remove correlations between the Breit-Wigner function and the binned 
S wave amplitude. The Kn P wave in Eq. (18) is represented by the JT*(892) and i^*(1680) 
Breit-Wigner functions, where the Cx*(i680) parameters are allowed to float in the fit. As 
usual cx*(892) = 1, and all other magnitudes and phases are defined with respect to fT*(892). 
The fT*(1410) is excluded as mentioned previously. The Kn D wave in Eq. (19) is repre- 
sented by the ^(1430) Breit-Wigner function with its parameters ck*(uso) allowed to float 
in the fit. The 1 = 2 7r + 7r + S wave is represented by the unitary amplitude from Eq. (26) 
fixing the m min and m max parameters to their optimal values from model 12. The a nn and 
parameters are allowed to float in the fit. 

The results of the QMIPWA fit are shown in Table VII for resonance parameters and 
Table VIII for the Kit S wave. Figure 10 shows the Dalitz plot projections of this fit. The 
measured Kn S wave and their comparison with model 12 components are displayed in 
Fig. 11. Both the magnitude and phase are different from those of model 12. 

Table VIII (right-hand side) also shows the result of a similar fit for a total binned S 
wave. For this fit the 1^(1430) resonance is accounted for in the binned S wave, and all 
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TABLE VII: Fit results for model 12, with a Breit-Wigner function or with a complex pole for the 
k, and QMIPWA. The FF are shown for a single Kir wave and need to be doubled as indicated by 
the "2x" symbol in row titles. 



Mode 


Parameter 


Model 12 (B-W for k) 


Model 12 


QMIPWA 


7T(892W+ 


a 


1 - fixed 


1 - fixed 


1 - fixed 




(°) 


- fixed 


- fixed 


- fixed 




FF (%) 2x 


5.15±0.24 


5.27±0.08±0.15 


4.94+0.23 




m (MeV/c 2 ) 


895.4±0.2 


895.7±0.2±0.3 


895.7 - fixed 




T (MeV/c 2 ) 


44.5±0.7 


45.3±0.5±0.6 


45.3 - fixed 


iT(1680)7r+ 


a 


4.45±0.23 


3.38±0.16±0.78 


2.88+0.84 




(°) 


43.3±3.6 


68.2±1.6±13 


113+14 




FF (%) 2x 


0.238±0.024 


0.144+0.013+0.12 


0.098+0.059 


^2(1430)7r+ 


a 


0.866±0.030 


0.915+0.025+0.04 


0.794+0.073 




4> (°) 


-17.4±3.5 


-17.4+2.3+2.0 


14.8+9.0 




FF (%) 2x 


0.124±0.011 


0.145±0.009±0.03 


0.102+0.020 


iTn(1430)7r+ 

u v / 


a 


3.97±0.15 


3.74+0.02+0.06 


3.74 - fixed 




<M°) 


45.1±0.9 


51.1+0.3+1.6 


51.1 - fixed 




FF (%) 2x 


7.53±0.65 


7.05+0.14+0.55 


6.65+0.31 




m (MeV/c 2 ) 


1461.1±1.0 


1466.6+0.7+3.4 


1466.6 - fixed 




T (MeV/c 2 ) 


177.9±3.1 


174.2+1.9+3.2 


174.2 - fixed 


KTT + 


a 


5.69±0.17 


10.80+0.05+0.35 







(°) 


-149.9±1.2 


148.4+0.3+1.6 







FF (%) 2x 


8.5±0.5 


21.6+0.3+3.2 





Pole 


5ftm (MeV/c 2 ) 
9m (MeV/c 2 ) 




706.0+1.8+22.8 
-319.4+2.2+20.2 




Breit-Wigner 


m (MeV/c 2 ) 
T (MeV/c 2 ) 


888.0±1.9 
550.4±11.8 






NR 


a 


17.1±0.4 


23.3+0.1+1.6 







(°) 


1.9±1.7 


29.7+0.2+3.0 







FF (%) 


38.0+1.9 


73.8+0.8+9.6 





Binned K~tt + S wave 


a 








See 




(°) 








Table VIII 




FF (%) 2x 








41.9±1.9 


1 = 2 7r + 7r + S wave 


a 


30.3+2.7 


25.5±0.3±2.9 


33.1±2.6 




(°) 


86.3+3.3 


75.4+0.6+10 


66.2+3.5 




FF (%) 


13.4+2.3 


9.8+0.2+2.0 


15.5+2.8 


Equation (28) 


At/ 


1 


1 


1 




m min (MeV/c 2 ) 


1265+8 


1256±5±4 


1256.3 - fixed 




TOmax (MeV/c 2 ) 


1529+31 


1498+13+30 


1498.2 - fixed 


Form factor 


^(GeV- 1 ) 


5 - fixed 


5 - fixed 


5 - fixed 




^at,e(GeV _1 ) 


- fixed 


- fixed 


- fixed 


Other R — ► Ktt 


^(GeV- 1 ) 


1.5 - fixed 


1.5 - fixed 


1.5 - fixed 




E FF, (%) 


94.4 


152.0 


122.8 


Goodness 


x 2 /^ 


426/385 


416/385 


359/347 




Probability (%) 


7.4% 


13.2% 


31.5% 
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other parameters are fixed to their values from the nominal QMIPWA fit. 

B. Cross-check for binned 1 = 2 7T+7T+ S wave 

In this approach we also check how much the 1 = 2 7r + 7r + S wave might be different 
from its analytic approximation in model 12. The analytic expression, given in Eq. (26), is 
replaced by the binned wave from Eq. (29). We fix all 20 fit parameters to their values in 
model 12. Then, we fix to zero the magnitude of the analytic 1 = 2 7t + tt + S wave, a n7T = 0, 
and add a binned 1 = 2 7r + 7r + S wave; the magnitude and phase in each of the 18 bins of 
mi + „ + are allowed to float. Results of this fit are shown in Table IX. The measured 1 = 2 
7T + 7r + S wave is compared to the model 12 analytic function in Fig. 12. The change of the 
X 2 in this fit compared to model 12, A% 2 = 390 — 416 = —26, does not show a significant 
improvement for the change Ais = —36 + 20 = —16 in degrees of freedom due to fixing 20 
of the original parameters and introducing 36 new parameters. The data do not prefer the 
binned amplitude to the analytic model for the 1 = 2 S wave. 

C. Cross-checks for binned Ktv P and D waves 

To check how much the P and D waves might be different from their model parametriza- 
tion we use the same binned technique. All of the fit parameters are fixed to their optimal 
values from model 12. The P wave binned amplitude substitutes for the smooth K*(1680) 
resonance only. The sharp shape of the K*(892) resonance is accounted for by a Breit- 
Wigner function. The D wave binned amplitude substitutes for the ^(1430) contribution, 
even though it is not particularly smooth. The magnitude and phase in 26 bins of m 2 KlT are 
allowed to float for the P or D wave in two separate fits, respectively. The resolution for the 
D wave is worse, and we use only 18 bins in the m 2 Kw range [0.9, 2.7] (GeV/c 2 ) 2 . Results of 
these fits are shown in Table X. The measured P and D binned waves and their comparison 
with the model 12 components are displayed in Figs. 13 and 14, respectively. 

The relative fraction of the P and especially the D wave is much smaller than the S wave. 
This explains the poor resolution for these binned amplitudes. We find that our binned P 
and D waves are consistent with their substituted components in model 12. This cross-check 
adds some confidence to this quasi-model-independent technique. 

VI. CROSS-CHECKS AND SYSTEMATIC UNCERTAINTIES 
A. Systematic uncertainties 

In order to estimate systematic uncertainties of the fit parameters, we apply numerous 
variations to the fitting procedure and examine the change of the fit parameters. Because of 
the numerous variations, a quadratic sum of the resulting small changes due to the variations 
would lead to a systematic error that is too conservative. Instead, we treat all of the resulting 
changes from the variations as a frequency distribution and calculate the mean and rms of 
the distribution. The total systematic uncertainty is obtained as a quadratic sum of the 
difference between the obtained mean and the default value of that parameter, and the rms. 
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TABLE VIII: QMIPWA: Results for Kir S wave. Figure 11 shows the binned Kit S wave without 
Kq(1430) by the dots with error bars and the total Kir S wave by the solid curve. 



Bin 

# 

7T 


< (GeV/c 2 ) 2 
Bin range 


Binned Kit S wave w/o i^ *(1430) 


Total Kit S 


wave 


Magnitude (a.u.) 


r^nase ( ) 


Magnitude (a.u.) 


Phase (°) 


i 


n /i 

0.4 


- 0.5 


on oo_i_n cn_i_n 
zU.2o±0.o0±0.0l 


T1 0_l_/1 0_l_ A 1 

— j l.o±4.y±4.1 


19.29±0.56 


-66.0±3.2 


2 


n r 

0.5 


- (J.O 


on nn i n no \ r\ r A 

20.90±0.68±0.54 


-61.1±4.6±z.3 


20 38±0 51 


-54.3±2.8 


3 


0.6 


— U.7 


20.58±0.63±0.b9 


A o o i i r\ i o o 

-48.8±4.0±3.z 


20 35±0 45 


-40 9±2 2 


4 


n ^7 

0.7 


- (J. 8 


on /^o i n c /i i n c; 1 

z0.bz±0.b4±0.51 


-47.U=bz.4±U.7 


20 52±0 59 


-39 7±0 9 


5 


0.8 


— U.y 


on nc i n ^7/1 i 1 oo 

z0.95±0.74±l.oz 


A A 1 1 9 1 9 C 

-44.1±o.b±o.o 


20.84±0.49 


-36. Oil. 4 


6 


n n 
0.9 


1 n 
— l.U 


i n n7J_n £;q_i_i 1 /I 

19.9/±0.bo±1.14 


— oo.4±4.U±D.4 


19 97±0 36 


-28.7±2.7 


7 


1.0 


- 1.1 


in oc_i_n coin nc 

19.3b±0.58±0.96 


o/? 1 1 ') r_|_/f o 


19 97±0 34 


-15.7±2.6 


8 


1.1 


- l.Z 


1 ^7 oi _i_n coin t\a 

17.81±0.53±0.9b 


-lo.D±O.U±D.l 


18 85±0 30 


-6.9±2.0 


9 


1 o 
1.2 


1 9 

- l.o 


1 ^7 7n i n /in_i_n nn 
1 (. f 0±0.49±0.90 


— lo. / ±z.o±o.l 


19 05±0 27 

-1 ■ ' • V / ".y — 1 — V / . AJ 1 


-1.2±1.6 


10 


l.o 


1 A 

- 1.4 


1 / z±0.4/±1.05 


-o.o±z.o±o.o 


19.47±0.25 


4.8±1.4 


11 


1.4 


- 1.5 


1*7 i Q_i_n /ir._i_i nc 
17.1o±0.4o±1.05 


9 0_LO A A- A T 
— O.Z±Z.4±4. / 


19.36±0.23 


11.6±1.3 


12 


1.5 


- l.b 


17.1b±0.45±1.04 


n r i o o i r /? 


19.80±0.22 


17.1±1.3 


13 


l.o 


-1.7 


17.09±0.4b±1.14 


o.4±z.o±4.U 


20.35±0.22 


24.3±1.3 


14 


1. / 


- l.o 


1 n/i_i_n /i^7_i_i n/^ 
lb.94±U.4/±l.Ub 


T FT 1 O Q_L/1 ^ 


20.61±0.24 


30.1±1.3 


15 


l.o 


i n 
- 1.9 


1 A 1 in /IOI1 nc 

lb.41±0.48±1.05 


n c i o /i i o t 

9.5±z.4±J.7 


20 38±0 27 


38 Oil 3 


16 


1.9 


- 2.0 


i c ni _i_n ci _i_n ^7C 
15.91±0.51±0. /5 


1 O Q_LO 7_LQ O 

lz.o±z. f ±O.Z 


19 86±0 34 


49.1±1.4 


17 


2.0 


- 2.1 


15.97±0.5b±l.ll 


16.9±2.9±3.1 


19 07±0 40 


63.4±1.4 


18 


2.1 


- 2.2 


15.72±0.61±1.11 


15.7±2.5±2.2 


14.20±0.40 


76.6±1.9 


19 


2.2 


- 2.3 


16.54±0.63±1.17 


17.3±2.1±2.1 


9.29±0.39 


75.2±2.6 


20 


2.3 


- 2.4 


16.64±0.71±1.42 


20.1±2.1±3.4 


7.10±0.43 


57.5±3.1 


21 


2.4 


- 2.5 


15.91±0.86±1.61 


19.8±2.3±2.0 


6.75±0.43 


33.5±3.2 


22 


2.5 


- 2.6 


17.25±1.03±1.83 


21.5±2.4±2.0 


9.46±0.51 


26.0±2.5 


23 


2.6 


- 2.7 


17.24±1.21±2.38 


24.1±2.6±2.9 


10.74±0.55 


26.7±2.6 


24 


2.7 


- 2.8 


17.59±1.34±2.28 


29.0±2.8±1.6 


12.04±0.64 


31.8±2.5 


25 


2.8 


- 2.9 


16.51±1.65±2.56 


32.7±3.5±3.5 


11.79±0.79 


36.0±3.3 


26 


2.9 


- 3.0 


14.07±3.28±3.32 


35.6±6.0±3.8 


9.75±2.50 


38.6±6.1 



We consider the six divisions of our sample: splitting events evenly between earlier and 
later datasets; D + and D~ decays; tight [la(AE) x la(mBc)]; an d loose [3a(AE) x3cr(m B c)] 
signal box cuts. Other systematic variations common to all of the models (model C, model 12, 
and QMIPWA), are as follows: 

• Float / S ig. We find that it is always consistent with its nominal value of 0.9892. 

• Drop the background term and set / S j g = 1. 

• Float the efficiency coefficients in a simultaneous fit to data and MC samples. In this 
case variations of the efficiency parameters are constrained by the MC sample. 

• Improve the precision of the calculation of the normalization integrals from Eq. (3) by 
an order of magnitude. 
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TABLE IX: QMIPWA: Results for 7r + 7r + S wave, which are also shown by dots with error bars 
in Fig. 12. Variation of the x 2 an d number of degrees of freedom, Au, is shown with respect to 
model 12. 



Bin 


(GeV/c 2 ) 2 


1 = 2 7T + 7T+ 


5 wave 


# 


Bin range 


Magnitude (a.u.) 


Phase (°) 


1 


0.1 — 0.2 


3.62±0.44 


-113.9±7.0 


2 


0.2 — 0.3 


4.31±0.37 


-115.8±4.7 


3 


0.3 — 0.4 


4.51±0.33 


-120.7±3.6 


4 


0.4 — 0.5 


6.17±0.33 


-116.8±2.4 


5 


0.5 — 0.6 


6.84±0.34 


-118.9±2.1 


6 


0.6 — 0.7 


8.13±0.35 


-121.8±1.8 


7 


0.7 — 0.8 


7.77±0.35 


-121.7±1.9 


8 


0.8 — 0.9 


8.65±0.35 


-125.5±1.7 


9 


0.9 — 1.0 


8.95±0.37 


-126.3±1.6 


10 


1.0 — 1.1 


9.61±0.42 


-127.9±1.5 


11 


1.1 — 1.2 


11.69±0.55 


-127.9±1.3 


12 


1.2 — 1.3 


10.04±0.93 


-132. Oil. 5 


13 


1.3 — 1.4 


12.43±1.01 


-131. Oil. 2 


14 


1.4 — 1.5 


12.92±0.98 


-131.1±1.3 


15 


1.5 — 1.6 


9.51±0.98 


-131.4±1.8 


16 


1.6 — 1.7 


11.56±1.00 


-134.7±1.8 


17 


1.7 — 1.8 


11.75±1.10 


-142.5±2.1 


18 


1.8 — 1.9 


10.52±2.02 


-156.1±5.2 




FF(%) 


9.8±0.4 




A x 2 /Av 


(390 - 416)/(- 


36 + 20) 



We repeated the fits with the radii for the Blatt-Weisskopf [19] form factors a factor of 
2 larger and smaller than their nominal values and found negligible change in the central 
result. 

Depending on the model we apply additional systematic variations, which are also in- 
cluded in systematic uncertainties. In model C and model 12, discussed in Sec. IV, we 
consider variation of the resonance parameters as follows: 

• Float parameters of the if* (892) and 1^(1430) resonances. The mass and width are 
shown in Table IV. 

• Float parameters of the iQ(1430) or if* (1680) resonances. The mass and width are 
shown in Table IV. 

• Add a if*(1410) contribution. We find that it is not significant: F Fk* (mo) < 0.1% at 
90% C.L. 

• In model C we test the event sample selected without a kinematic fit with appropriate 
coefficients for efficiency and background shapes. We find that the resulting variation 
is negligible, and this study is not performed for the other fit models. 



23 



TABLE X: Results for Kit P and D binned waves using model 12, which is also shown by dots with 
error bars in Figs. 13 and 14, respectively. Variation of the \ 2 an d number of degrees of freedom, 
Au, is shown with respect to model 12. The D wave is used only in the range [0.9,2.7] (GeV/c 2 ) 2 
as explained in the text. 



Bin 


(GeV/c 2 ) 2 


Binned P wave for IT" (1680) 


Binned D wave for K|(1430) 


# 


Bin range 


iviagnuuQe ^a.u.j 


r nase ^ ) 


iviagmcuoe i^a.u.j 


c nase ^ ) 


1 


n a 

U.4 — 


U.o 


1 nft-H OR 

i.yo±i.zo 


a i n-i-oa n 

41.U±Zo.U 






2 


n k 

U.O — 


n ft 
U.O 


z.yu±i.iz 


1 no 0-L1 1 o 
iyy.z±n.z 






3 


n ft 

U.O — 


n 7 
U. 1 


n ft7_i_n hq 
U.O ( ±u.yo 


ou.i±io.o 






4 


n 7 

u. / — 


n 8 
U.o 


n 81 xn 
U.olzhU.oO 


Ad \ I 1 8 /I 
40.4±lo.4 






5 


u.o — 


n q 
u.y 


U.OUztU.oO 


lUZ.oztoo. / 






6 


n o 

U.y — 


1 n 
l.U 


n c: Q_|_ 1 oft 

U.Oozhl.oO 


7*3 O-I-Of; 7 
/ o.ZzhZO. I 


n 71 -i~n 1 7 
U. i l±u.l ( 


ft1 S-I-OQ ft 

— 01.o±zo.0 


7 


i n 

l.U — 


1.1 


1 Q7_l_n 77 

l.o i ±U. ( ( 


7n ^-L7 m 

<y.o± i .ui 


U.ol±U.l4 


k nj-ftn ft 

— o.u±oy.o 


8 


1.1 — 


L.Z 


n nnxn /io 
U.Uy±U.4z 


K1 7-LOO /I 

01. / ±zz .4 


n ftft_i_n i ft 
U.OOztU.lO 


1 K ft_l_07 O 
— 10.0±Z 1 .Z 


9 


1 

l.Z — 


l.o 


i i o_i_n q c 
l.lz±U.oo 


81 q_i_i n ft 
ol.oztlU.O 


ii/i _i_n i n 
1.14±U.iy 


Kd n_i_ ion 
— OU.Uztlz.y 


10 


l.o — 


1 /I 
1.4 


1 /IO-l-n QQ 

1.4y±U. oo 


7/1 A 1 8 70 

/ 4.4±o. ( y 


i nn-i~n i o 

i.uy±u.iy 


08 74-1 Q 1 
— Zo. ( ±lo.l 


11 


1 /I 

1.4 


1.0 


n 81 4-0 

U.olzhU.oo 


8Q O-l-l Q Q 

oo.y±io.o 


n o74~n i o 

u.y / ±u.iy 


n n_i_i /I o 

— y.y±i4.z 


12 


l.o — 


1 £ 
1.0 


l.Z0±U.oo 


on 8-L1 n ft 
oU.o±lU.O 


i onxn oi 
l.oUztU.Zl 


1 ^ ^j-1 A n 
— 10.0±14.U 


13 


1 £ 

l.o — 


1 7 

1. / 


i Qn_i_n Q/i 
l.oU±U.o4 


7Q Q_l_n 70 

1 o.o±y. < y 


i A £_i_n o/i 
1.40±U.z4 


11 ft_l_ 1 Q K 

ll.b±lo.O 


14 


i 7 


1 8 

-L .O 


1 77_l_n Qft 

L. 1 1 ±U.o0 


77 7_1_8 Q/1 

M . ( ±o.o4 


o i c_i_n oft 
z.loztU.ZO 


1ft n o 
— 10.o±lU.z 


15 


1.8 - 


1.9 


i i ft_i_n An 


QO 111 1 1 
OO . 1 _L_ J. _L . _L 


q oe_|_n oo 


94 44-7 


16 


1.9 - 


2.0 


1.79±0.45 


90.9±9.91 


5.77±0.38 


31.9±5.45 


17 


2.0 - 


2.1 


1.86±0.42 


102.1±11.2 


4.77±0.52 


70.2±5.90 


18 


2.1 - 


2.2 


2.42±0.45 


104.3±12.2 


4.45±0.70 


100.2±7.41 


19 


2.2 - 


2.3 


2.58±0.52 


115.7±12.6 


4.69±0.97 


127.8±8.52 


20 


2.3 - 


2.4 


1.66±0.55 


119.6±16.6 


3.79±1.13 


116.8±16.3 


21 


2.4 - 


2.5 


3.76±0.60 


108.8±11.4 


0.18±1.42 


-3.3±84.7 


22 


2.5 - 


2.6 


4.19±0.75 


113.6±13.7 


4.86±1.72 


135.8±33.2 


23 


2.6 - 


2.7 


6.71±0.98 


105.0±16.9 


1.54±2.58 


184.8±78.5 


24 


2.7 - 


2.8 


3.35±1.59 


142.4±29.4 






25 


2.8 - 


2.9 


7.03±2.49 


183.2±22.5 






26 


2.9 - 


3.0 


32.66±17.9 


232.2±11.6 








FF(%), 2x 


0.20±0.03 


0.15±0.01 




A x 2 /Ais 


(373 -416)/(- 


■52 + 20) 


(400 - 416)/(- 


-36 + 20) 



• In QMIPWA we use the binned K tt + S wave as a histogram without linear interpo- 
lation. 

None of these variations reveal an obviously dominant source of systematic uncertainty. 



B. Additional cross-checks 



Common cross-checks listed below bring us to models that are different enough that we 
do not include these results in the total systematic uncertainty. 

o Remove the if £ (1430) resonance contribution. We find a significant degradation of 
the fit quality. For example, A% 2 = 370 in model C. 

o Remove the if* (1680) resonance contribution. We find a significant degradation of 
the fit quality. For example, A% 2 = 250 in model C. 

Based on these fits we conclude that the ^(1430) and if* (1680) resonances cannot be 
removed. Their contributions are small (FF< 1%) but significant and the fit quality is 
very poor without these resonances. We apply a similar backward check for model 12 and 
QMIPWA: 

o We remove the 1 = 2 n + n + S wave contribution and get a poor quality fit. 
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FIG. 8: Projections of the fit to the Dalitz plot with model C on (a) m 2 (KTr) (two entries per 
event) and (b) m 2 (7T7r) variables. Residuals between data and the total p.d.f. are shown by dots 
with statistical error bars on the top insets along with minor contributions from the if* (1680) and 
if 2 (1430) resonances, plotted with factor x4. 
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FIG. 9: Projections of the fit to the Dalitz plot with model 12 on (a) m 2 (K-ir) (two entries per 
event) and (b) m 2 (7T7r) variables. Residuals between data and the total p.d.f. are shown by dots 
with statistical error bars on the top insets along with minor contributions from the If* (1680) and 
if* (1430) resonances, plotted with factor xlO. 



VII. DISCUSSION 

A. if* (1430) and if* (892) parameters 

The ifg(1430) parameters, which we obtain in fits with the isobar models, listed in 
Table IV as well as in Tables V-VII, are consistent with each other and are significantly 
different from the PDG [1] values which are dominated by the LASS measurement. Our 
data prefer a Aq(1430) resonance that is 50 MeV/c 2 heavier and about 2 times narrower. 
Our result is consistent with E791 [6] and FOCUS [23] measurements. Possible explanations 
include that the Aq(1430) resonance parameters depend on its production mechanism which 
are different in Kir scattering and D-meson decays or that the Aq(1430) parameters are 
strongly model-dependent. In particular, they may depend on interference with other S 
wave contributions. 

In all of our fits, as seen in Tables IV and VII, the A"* (892) parameters are consistent with 
each other. The measured masses are consistent with the nominal value [1], but the widths 
are about 5 MeV/c 2 narrower. The A*(892) width in the PDG is an average value over 
about 20 experimental measurements. Our measurement does not contradict any single 
result included in PDG. The largest difference is about 3 standard deviations from the 
LASS value. We also tested the A* (892) width model dependence. The largest variation is 
expected from applying Blatt-Weisskopf form factors. A variation of the radial parameters 
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FIG. 10: Projections of the fit to the Dalitz plot with QMIPWA on (a) m 2 (K7r) (two entries per 
event) and (b) to 2 (7T7t) variables. Residuals between data and the total p.d.f. are shown by dots 
with statistical error bars on the top inset along with minor contributions from the X*(1680) and 
if|(1430) resonances, plotted with factor xlO. 

from zero to 4 times their nominal values causes a variation in the width of 0.7 MeV/c 2 , which 
cannot explain the large difference with the PDG average. Recent results from FOCUS [24] 
and Belle [25] also indicate a if* (892) width smaller than the PDG value. 



B. Partial waves 

The factorized Gaussian form factors used in model C for the Kir S wave components may 
cause an enhancement of the complex function magnitude at low Kir mass, as mentioned in 
Sec. IV B. However, the form factor is a real function and does not change the phase of the 
complex S wave. In model 12 and QMIPWA we measure the total Kir S wave amplitude 
without using a form factor. It means that the measured S wave absorbs the form factor. 
The measured S wave magnitude is essentially constant up to 1.4 GeV/c 2 , as demonstrated 
in Fig. 11. The phase of the complex S wave amplitude shows smooth variation from -80° 
at Kit threshold to 40° at 1.4 GeV/c 2 . At higher mass the amplitude is distorted by a 
contribution from the Kq(1430) resonance. Our data require a dominant contribution from 
the Kn S wave, which at low mass is not well described as a regular resonance structure. 

In QMIPWA the S wave is measured in a model-independent way, while the P and D 
waves are parametrized by Breit-Wigner resonances. In cross-checks we add more freedom 
to the P and D waves by replacing the if* (1680) and ifj;(1430), respectively, Breit-Wigner 
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FIG. 11: The magnitude and phase of the Kir S wave in model 12 and QMIPWA. The dots with 
error bars for statistical uncertainties and the solid curve show the binned Kir S wave without 
ifg(1430) and the total Kir S wave from Table VIII, respectively. Other curves show the S wave 
components of model 12 with parameters from Table VII. 

shapes with the binned amplitudes, as illustrated in Figs. 13 and 14. We did not find 
improvement in the fit quality or any significant deviation of the binned wave from the 
analytic function of the isobar model. The systematic uncertainties described above are 
larger than variations caused by the possible model dependence of the P and D waves. 

In model 12 and QMIPWA we find that the 1 = 2 ir + ir + S wave amplitude is not consistent 
with a constant term as seen in Fig. 12. Its behavior is well modeled by the analytic function 
[Eq. (26)]; the binned wave also describes the behavior consistently. 

C. Comparison of fit models 

Note that in model C, model 12, and QMIPWA (Tables IV- VII) the magnitudes and 
phases for the resonance contributions if* (892) (by construction), if,* (1430), iff (1430), and 
if* (1680) are consistent. On the other hand, the fit fractions for the resonances differ 
significantly among the models. The differences in fit fractions arise from our models of 
the Kir S wave and the 1 = 2 ir + ir + S wave when it is included. Indeed, in model C the 
Kir S wave amplitude is represented by the NR term and a broad Breit-Wigner resonance 
for the k. They have small magnitudes, and they interfere constructively. This leads to 
small fit fractions and the sum of all fit fractions is less than 100%. In model 12 the Kir S 
wave amplitude is represented by the NR term and the complex pole for the k. The i = 2 
ir + ir + S wave is also included in this model. These amplitudes are large, and they interfere 
destructively. This leads to large fit fractions, and the sum of all fit fractions is greater than 
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FIG. 12: The magnitude and phase of the 1 = 2 tt + tt + S wave in model 12 and QMIPWA. The dots 
with error bars for statistical uncertainties represent results from Table IX. Other curves show the 
1 = 2 7T + 7r + S wave of the model 12 and QMIPWA with parameters from Table VII. 



100%. Both models describe the total Kir S wave with complex functions that qualitatively 
show very similar behavior. However, the fit fractions, especially for the components of the 
Kn S wave, strongly depend on the assumed composition of this amplitude. 

In contrast the QMIPWA represents the Kn S wave amplitude as a single function. We 
make no assumptions about the composition of the Kn S wave in this case. It also includes 
the 1=2 7r + 7r + S wave, which plays a key role in improving the fit quality. We feel this 
approach is the most reliable presentation of results for the fit fractions for resonance and 
S wave contributions to the D + — ■> K~it + -n + decay. 



VIII. SUMMARY 

We describe a partial wave analysis of the D + — > K~ir + ir + events on the Dalitz plot. 
We use the CLEO-c data set of 572 pb^ 1 of e + e~ collisions accumulated at the i/>(3770), 
which corresponds to a sample of 1.6 x 10 6 D + D~ pairs produced in the process e + e _ — > 
^(3770) — > D + D~ . We select 140793 candidate events for the Dalitz plot with a small 
background of 1.1%. We compare our results with the best previous measurements from 
E791 [6] using the isobar model. Our results agree with the E791 measurement, as shown 
in Table V for their best model. The fit quality can be improved if we add the 1 = 2 -k + ti + 
S wave contribution, as presented in Table VII. 

We apply a model-independent approach, developed in Ref . [7] , to measure the magnitude 
and phase of the Kn S wave in the invariant mass range from threshold to the maximum 
value in this decay. In contrast to E791, we have measured the Kn S wave without factor- 
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FIG. 13: The magnitude and phase of the Kir P wave in model 12 and QMIPWA. The dots with 
error bars for statistical uncertainties represent results from Table X. Other curves show the P 
wave components of model 12 and QMIPWA with parameters from Table VII. 



ization of the form factor for scalar resonances. Our results on the Ktt S wave phase and 
magnitude measurement are presented in Table VIII and in Fig. 11. The accuracy of the 
Ktt S wave measurement is improved, compared to the only previous measurement [7] . We 
find that the total observed S wave magnitude in the D + — > K~tt + tt + decay is essentially 
constant from Ktt production threshold to 1.4 GeV/c 2 . The phase shows smooth variation 
from -80° to 40° in the same range. At higher invariant mass m Kn > 1.4 GeV/c 2 , the S wave 
behavior is dominated by the Kg (1430) resonance. We find that the P wave contribution is 
dominated by K*(892) and X*(1680) Breit-Wigner resonances, and the D wave has only a 
contribution from i^2(1430). Using binning techniques, we find no significant deviation of 
the P and D waves from the isobar model, as demonstrated in Figs. 13 and 14. 

In the model-independent approach for the 1 = 2 it + it + S wave we obtain the binned 
amplitude parameters, listed in Table IX and shown in Fig. 12, which are consistent with the 
analytic form of this wave. We find that the 1 = 2 it + it + S wave has a nonuniform variation 
in the amplitude across the m(7r + 7r + ) kinematic range with a fit fraction of 10%-15%. As 
expected, the measured amplitude behavior and the fit fraction of the Ktt S wave changes 
slightly with the addition of the 1 = 2 tt + tt + S wave. The addition of the 1 = 2 it + it + S 
wave to either the isobar model or the model-independent partial wave approach is the key 
piece that gives good agreement with the data in both cases. 
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FIG. 14: The magnitude and phase of the Kir D wave in model 12 and QMIPWA. The dots with 
error bars for statistical uncertainties represent results from Table X. Other curves show the D 
wave components of model 12 and QMIPWA with parameters from Table VII. 
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APPENDIX A: APPENDIX: KINEMATIC VARIABLES AND ANGULAR DIS- 
TRIBUTIONS 

Following previous CLEO analyses (see, for example, Ref. [18]) we use the angular dis- 
tributions [Eqs. (10)— (12)] obtained from a covariant-tensor formalism. The E791 form 
[7] was applied to the orbital momentum partial wave decomposition using the term 
(-2P a P c ) L P L (cosf^) in stead of il L (s,t) in Eq. (9). Here V L (cos 6) is a Legendre polynomial, 
P a , P c , and 9 are the momenta and the angle between particles a and c in the resonance 
R rest frame, respectively. In this section we show that both approaches are equivalent up 
to constant coefficients. As a by-product we also obtain expressions for the decay momenta 
used in the Blatt-Weisskopf form factors [Eqs. (13)-(16)]. 

The kinematic variables in the decay under study, schematically shown in Fig. 6, can 
be expressed in terms of invariant variables (masses and invariant masses) in the decaying 
particle d (D + meson) or resonance R rest frames. The energy and momentum of particle a in 
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the resonance R or (ab) rest frame can be obtained from the 4-momentum balance equation 
Pb = PR~Pa, where p R = (m ab , 0) and p a = (E a , P a ). Then, p 2 b = m 2 b = m 2 ab + m 2 a - 2m ab E a , 
giving the energy 

E a = <^<~ m2 \ (Al) 
2m ab ' 1 } 

and relevant momentum squared 

?l = El-ml = \( ml - 2ml - 2m\ + ^Z^t) . (A2) 



m 



ab 



The energy and momentum of the particle c in the resonance R rest frame can be obtained 
from the particle d invariant mass squared, m 2 d = (p c + Pr) 2 . In the resonance rest frame 
p c = (E c , P c ) and Pr = (m ab , 0), and therefore m 2 d = m 2 + m 2 ab + 2m ab E c , from which we 
find its energy 



ml — m nh — m f . 
E = _A ± c A3 

2m ab ' 1 ; 



and the associated momentum squared 



P, 2 = ti -m 2 c = \ (m 2 ab - 2m 2 d - 2m\ + K ) ( A I ) 

1 ~ rn ab 



The energy and momentum of the same particle c in the D meson rest frame, denoted 
here by the asterisk, have different expressions: 

= m 2 -m 2 ab + m 2 2 = 2 _ ^ 

2m d 

The angular distributions in the E791 [7] analysis are defined by the Legendre polynomials 

V L (cos9), 

Po(x) = l, V 1 (x)=x, V 2 (x) = ^—, V 3 (x) = (A6) 

where 9 is the angle between particles a and c in the resonance R rest frame. This cos 6* 
can be expressed through the known energies and momenta of particles a and c and their 
measured invariant mass squared m 2 c = (p a + p c ) 2 = m 2 + m 2 + 2E a E c — 2P a P c cos 9, so 

cos9= < + m c + ^ Ec - m - . (A7) 
2P a P c 1 ' 

Substituting E a and E c from Eqs. (Al) and (A3) in the numerator of Eq. (A7) we get 
mq/3 _ 1 ( m 2 m 2 | (m 2 d -m 2 c )(m 2 a -m 2 ) \ 

** r a r c V "''ab 7 

Note that these angular distributions are equivalent to formulas used in Eqs. (10)— (12) up 
to constant factors. Indeed, comparing the expressions for (— 2P a P c ) L Pi(cos6 l ) from Ref. [7] 
with flj,we get 
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(-2P a P c )°P o (cos0) 
(-2P a P c ) 1 P 1 (cos^) 



-2P a P c cos# 



1 




(A9) 
(A10) 



(-2P a P c ) 2 P 2 (cos#) 



4(P a P c ) 



2 



(AH) 



[1] W.-M. Yao et al, Journal of Physics G 33, 1 (2006). 

[2] J. Adler et al. (MARK III Collaboration), Phys. Lett. B 196, 107 (1987). 

[3] M.P. Alvarez et al. (NA14/2 Collaboration), Z. Phys. C 50, 11 (1991). 

[4] J.C. Anjos et al. (E691 Collaboration), Phys. Rev. D 48, 56 (1993). 

[5] PL. Frabetti et al. (E687 Collaboration), Phys. Lett. B 331, 217 (1994). 

[6] E.M. Aitala et al. (E791 Collaboration), Phys. Rev. Lett. 89, 121801 (2002). 

[7] E.M. Aitala et al. (E791 Collaboration), Phys. Rev. D 73, 032004 (2006). 

[8] R.H. Dalitz, Philos. Mag. 44, 1068 (1953). 

[9] D. Bugg, hep-ex/0510014. 
[10] D. Bugg, Phys. Lett. B 632, 471 (2006). 
[11] J. A. Oiler, Phys. Rev. D 71, 054030 (2005). 

[12] W. Hoogland et al, Nucl. Phys. B 126, 109 (1977); N.B. Durusoy et al, Phys. Lett. 45B, 
517 (1973); B.S. Zou et al, published in Hadron 2003 Proceedings, AIP Conf. Proc. 717:347 
(2004) [hep-ph/0405118]. 

[13] G. Viehhauser, Nucl. Instrum. Methods A 462, 146 (2001); D. Peterson et al, Nucl. Instrum. 
Methods Phys. Res., Sect. A 478, 142 (2002); Y. Kubota et al, Nucl. Instrum. Methods Phys. 
Res., Sect. A 320, 66 (1992); R.A. Briere et al. (CESR-c and CLEO-c Taskforces, CLEO-c 
Collaboration), Cornell University, LEPP Report No. CLNS 01/1742 (2001) (unpublished). 

[14] S. Dobbs et al. (CLEO Collaboration), Phys. Rev. D 76, 112001 (2007). 

[15] H. Albrecht et al, Phys. Lett. B 229, 304 (1989). 

[16] G. Bonvicini et al. (CLEO Collaboration), Phys. Rev. D 76, 012001 (2007). 

[17] D.J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. A 462, 152 (2001). 

[18] S. Kopp et al. (CLEO Collaboration), Phys. Rev. D 63, 092001 (2001). 

[19] J.M. Blatt and V.F. Weisskopf, Theoretical Nuclear Physics, Wiley, New York, 1951, p. 361. 

[20] N.A. Tornqvist, Z. Phys. C 68, 647 (1995). 

[21] N.N. Achasov and G.N. Shestakov, Phys. Rev. D 67, 114018 (2003). 
[22] D.E. Groom et al, The European Physical Journal C 15, 1 (2000). 
[23] J.M. Link et al. (FOCUS Collaboration), Phys. Lett. B 653, 1 (2007). 
[24] J.M. Link et al. (FOCUS Collaboration), Phys. Lett. B 621, 72 (2005). 
[25] D. Epifanov et al. (Belle Collaboration), Phys. Lett. B 654, 65 (2007). 



33 



